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ABSTRACT 

We present a detailed spectral analysis of the data obtained from NGC 3783 during the period 2000-2001 using 
Chandra. The data were split in various ways to look for time- and luminosity-dependent spectral variations. 
This analysis, along with the measured equivalent widths of a large number of X-ray lines and photoionization 
calculations, lead us to the following results and conclusions. 1) NGC 3783 fluctuated in luminosity by a 
factor ^ L5 during individual observations (most of which were of 170 ks duration). These fluctuations were 
not associated with significant spectral variations. 2) On a longer time scale (20-120 days), we found the 
source to exhibit two very different spectral shapes. The main difference between these can be well-described 
by the appearance (in the "high state") and disappearance (in the "low state") of a spectral component that 
dominates the underlying continuum at the longest wavelengths. Contrary to the case in other objects, the 
spectral variations are not related to the brightening or the fading of the continuum at short wavelengths in any 
simple way. NGC 3783 seems to be the first AGN to show this unusual behavior. 3) The appearance of the 
soft continuum component is consistent with being the only spectral variation, and there is no need to invoke 
changes in the opacity of the absorbers lying along the line of sight. Indeed, we find all the absorption lines 
which can be reliably measured have the same equivalent widths (within the observational uncertainties) during 
high- and low-states. 4) Photoionization modeling indicates that a combination of three ionized absorbers, each 
split into two kinematic components, can explain the strengths of almost all the absorption lines and bound- 
free edges. These three components span a large range of ionization, and have total column of about 4 x 10^^ 
cm^^. Moreover, all three components are thermally stable and seem to have the same gas pressure. Thus 
all three may co-exist in the same volume of space. This is the first detection of such a multi-component, 
equilibrium gas in an AGN. 5) The only real discrepancy between our model and the observations concerns the 
range of wavelengths absorbed by the iron M-shell UTA feature. This most likely arises as the result of our 
underestimation of the poorly-known dielectronic recombination rates appropriate for these ions. We also note 
a small discrepancy in the calculated column density of O VI and discuss its possible origin. 6) The lower limit 
on the distance of the absorbing gas in NGC 3783 is between 0.2 and 3.2 pc, depending on the component of 
ionized gas considered. The assumption of pressure equilibrium imposes an upper limit of about 25 pc on the 
distance of the least-ionized component from the central source. 

Subject headings: galaxies: active — galaxies: individual (NGC 3783) — galaxies: nuclei — galaxies: Seyfert 
— techniques: spectroscopic — X-rays: galaxies 
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1. INTRODUCTION 

The barred-spiral galaxy NGC 3783 (V ~ 13.5 mag., z — 
0.0097) hosts a well-studied, type-I active galactic nucleus 
(AGN) with prominent broad emission lines and strong X- 
ray absorption features. The object has been observed ex- 
tensively with almost all X-ray instruments, most recently by 
Chandra (Kaspi et al. 2002, hereafter Paper I) and XMM- 
Newton (Blustin et al 2002). The 2-10 keV flux of NGC 3783 
varies in the range ^ (4-9) x 10^'^ ergscm^^s^', and its 
mean 2-10 keV luminosity is ^ 1 .5 x 10^^ erg s^' (for Hq = 
70 kms^' Mpc^' and qo = 0.5). Paper I gives an extensive 
list of references and a comprehensive summary of historical 
observations, including ground-based and UV (HST) observa- 
tions. It also discusses the unique Chandra data set obtained 
in 2000-2001 . These observations consist of a relatively short 
observation performed in 2000 January, and five longer ob- 
servations performed in 2001 February - June, separated by 
various intervals from 2 to 120 days (see §2). Paper I con- 
tains numerous illustrations of the mean spectrum, absorption 
line profiles and detailed measurements of many absorption 
and emission lines. Two other papers (Gabel et al. 2003a,b) 
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discuss the HST and FUSE data obtained as part of this muhi- 
waveband campaign. 

There have been several previous attempts to model the 
characteristics of the X-ray absorbing gas along the line of 
sight to NGC 3738. Here we comment only on the more de- 
tailed works. A very detailed work Kaspi et al (2001) at- 
tempted to fit the spectrum obtained with Chandra in 2000. 
The Kaspi et al. most successful model consists of two ab- 
sorbing shells, both outflowing from the central source with 
a velocity of ^ 600 km s~^ The gas in both shells was as- 
sumed to be turbulent with internal turbulent motion of ~ 250 
km s^' . The two shells had similar column densities, but dif- 
ferent ionization parameters. The model is consistent with 
the intensity and equivalent width (EW) of many (but not all) 
emission and absorption Unes observed in the spectrum. A 
major limitation of the Kaspi et al. (2001) work was the Um- 
ited signal-to-noise (S/N) of the spectrum, which results in 
large uncertainties on the model parameters. In addition, the 
limited duration of the only observation then available did not 
allow a meaningful analysis of any time-dependence. A new 
paper, by Krongold et al. (2003, hereafter K03), discussing 
the full Chandra data set, was accepted for publication after 
the submission of our paper. Some results of that work are 
relevant to our study and are discussed in the various sections 
below. 

This paper discusses the complete spectroscopic Chandra 
data of NGC 3783. We present our modeling of this material, 
with emphasis on the properties of the absorber(s) along the 
Une of sight. In §2 we describe the two X-ray spectral states 
discovered using the new observations. In §3 we explain the 
various ways we measured and modeled the spectrum. §4 con- 
tains a discussion of the new findings, and §5 summarizes our 
new results. 

2. THE TWO-STATE X-RAY SPECTRUM OF NGC 3783 

A full description of the Chandra data considered here is 
given in Paper I. In brief, there were 5 observations (each of 
~170 ks duration) obtained over the period 2001 February - 
June, and a shorter observation (of ^56 ks duration) obtained 
in 2000 January (and also discussed in Kaspi et al. 2001). All 
six observations were performed with the High-Energy Trans- 
mission Grating Spectrometer (HETGS) in place, and consist 
of a total exposure of 888.7 ks. All measurements in Paper 
I refer to a time-averaged spectrum, produced by combining 
the first-order spectra from both grating arms, using 0.01 A 
wide bins. Third-order data from the medium-energy grat- 
ing (MEG) were also used to compare the profiles of several 
short-wavelength lines with the first-order profiles obtained 
for several lines at longer wavelength. It was found that most 
of the resolved absorption lines (e.g. O VII, Ne IX and Si XIV 
Unes) consist of at least two kinematic components, outflow- 
ing at -500 ± 100 km s'^ and - 1000 ± 200 km s'^ (respec- 
tively). The overall absorption profiles in NGC 3783 covers 
the velocity range of to -1600 km s~^ (e.g. see Figs. 5 & 10 
in Paper I). 

A major goal of the present work is to investigate subsets 
of the Chandra data. We subdivided the data in various ways 
to search for spectral variations as a function of source lumi- 
nosity and/or time. This was done for both short (i.e. within 
the individual 170 ks observations) as well as for long (be- 
tween observations i.e. 20-90 days) time scales. The short 
time-scale behavior of NGC 3783 is illustrated in Fig.l which 
shows the short wavelength (2-10 A) flux as a function of time 
for aU observations biimed in intervals of 3170 sec in the left 



panels. The right panels shows the "softness" ratio, defined 
as the flux in the 15-25A band divided by that in the 2-10 A 
band, to illustrate now the different variability characteristics 
exhibited by the source above and below ~1 keV. 

From Fig. 1 it can be seen that the 2-10 A flux varied 
by about 50% over the 170 ks duration of each of the ob- 
servations. However, inspection of the corresponding soft- 
ness Ught curves shows that the spectral variations are much 
smaller and almost insignificant. Thus, the data show very 
little cahnges of the spectral energy distribution (SED) over 
time scales <4 days ("obsids" 2090 and 2091 are consecutive 
with total duration of about 4 days). To further test this find- 
ing, we have also extracted various spectra when the source 
was in different intensity states, but did not find any signifi- 
cant spectral variations during this period. We also note that 
Behar et al. (2003) did not find any spectral variations dur- 
ing a three-day observation of NGC 3783 in 2001 December 
using XMM-Newton. 

Contrary to the above, there is a very significant change of 
SED on longer time scales. In particular, there is a major 
change in the softness ratio between obsids 2092 and 2093 
(separated by 20 days) and obsids 2093 and 2094 (separated 
by 90 days). As shown in Fig. 1, obsid 2093 shows a higher 
mean count rate but exhibits a much softer spectrum. Direct 
comparison of the time-averaged spectra obtained during ob- 
sid 2093 with that obtained during obsid 2094 reveals the low 
energy part (A > 15 A) decreased in flux by a factor of ~ 4, 
but that the shorter wavelength continuum exhibited a much 
smaller decrease (factor ~ 1 .5) The opposite change occured 
between obsids 2092 and 2093. Most importantly, the soft- 
ness ratio variations are not simple luminosity-related effects. 
For example, the 2-10 A count rate in the middle of obsid 
2093 (Fig. 1) is a little lower than the 2-10 A flux toward 
the end of obsid 2094, yet the softness ratios are significantly 
different. A similar effect is seen in obsid 2093, during which 
the 2-10 A count rate varied by about 50% yet the softness 
ratio remaining approximately constant. 

In light of these results, we have divided the entire data set 
into groups with high and low softness ratios. We find four 
observations (2090, 2091, 2092 and 2094) with low softness 
ratio and two (0373 and 2093) with high softness ratio. Here- 
after we refer to these as the "low-" and "high-state" observa- 
tions, respectively. We find no significant differences between 
the mean spectra of the four individual low-state observations 
(except for a small intensity variations). Similarly we find no 
significant differences between the mean spectral shapes of 
the two high-state observations. 

Fig. 2 is a more detailed example of this phenomenon. It 
shows the softness ratio as a function of the 2-10 A flux for 
all the data using 3170s bins. A separation into two groups 
is apparent, and a standard KS statistical test confirms its sig- 
nificance. We have also investigated possible Unear correla- 
tions of the softness ratio with the 2-10 A flux within each 
of the states. For the low-state observations, we find a sig- 
nificant (> 99%) linear correlation between softness ratio and 
flux. However, as can be seen from Fig. 2, clearly this line 
does not connect the low- and high-state data. No significant 
correlation of softness ratio with flux was found for the high- 
state data. We conclude that the X-ray spectrum of NGC 3783 
fluctuates between two states of different softness ratio. The 
combined spectra of the two are shown in Fig. 3 and much of 
the rest of this paper is devoted to the analysis of this unusual 
behavior of NGC 3783. 
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Fig. 1.— MEG 2-10 A flux light curves (left) and softness ratios (ratio of the 2-10 and the 15-25 A count rates, right) for all Chandra observations of 
NGC 3783. The dashed lines across observations no. 2093 and 2094 represent constant flux (4.8) and softness ratio (0.15). They illustrate the fact that similar 
2-10 A fluxes (middle part of obs. 2093 and last part of obs. 2094) can be associated with very different softness ratios. 
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2-10 A Flux [10"" ergs cm"" s"' A"'] 

Fig. 2. — Softness ratio vs. 2-10 A flux for the data in Fig. 1. Tlie low 
and high-state observations are shown as filled and open symbols, respec- 
tively, and the various observations are denoted with different colors. The 
solid straight line is the linear regression fit to the low-state observations. No 
significant correlation of softness ratio vs. flux was found for the high state 
observations. 




wavelength (A) 



Fig. 3. — High-state (red, observations 0373 and 2093 combined) and low- 
state (blue, all other observations combined) spectra of NGC 3783 binned to 
0.04 A. The sohd fines are the two chosen intrinsic continua (see §3.2.5). The 
lower continuum is a single power-law with F = 1.6. The high continuum is 
made out of two power-laws with slopes F = 1.6 and F = 3.2 and relative 
normalization of 1:0.4 at 1 keV. These power-laws are discussed in §3.2.5. 



3. SPECTRAL ANALYSIS AND MODELING 
3. L Spectral differences between the high and the low states 

We have investigated the unusual variations observed in the 
spectrum of NGC 3783 in an attempt to understand their na- 
ture and their origin. In particular, we have attempted to an- 
swer the question of whether they are due to the response of 
the absorbing gas to the variations observed in the continuum. 

There are two ways to answer the above question - by a 
direct and detailed spectral comparison, and by modeling the 
two spectra trying to establish the origin of the differences. 
Fig. 4 is an example of the first approach. The diagrams 
shows a comparison of the low- and high-state spectra after 



Table 1. Equivalent Widths^ 



Ion & Line 


low state 


high state Kaspi et al. (2002) 


Si XIV A6. 182 


20.4 ±1.2 


22.6±1.7 


20.5 ±0.8 


Si XIII A6.648 


16.0±1.3 


14.6±1.3 


14.9 ±0.7 


Si XII A6.718'' 


2.9 ±0.5 


3.8±0.7 


3.0 ±0.4 


Si XI A6.778 


4.8 ±0.8 


5.9±1.2 


4.8 ±0.6 


SIX A6.859 


10.6±1.1 


9.4 ±1.4 


10.9±0.8 


SilX A6.931 


7.2±1.0 


7.7±1.6 


7.1 ±0.7 


Si VIII A6.999 


4.3 ±1.0 


4.4 ±1.4 


4.3 ±0.8 


S XVI A4.729 


11.5±1.6 


11.6±2.5 


10.7 ±1.2 


S XV A5.039 


8.7±1.7 


9.1±2.5 


9.2±1.2 



"Equivalent widths were measured as described in Paper I and are given in 
mA. 

''This line is blended with Mg XII A6. 738 with 1:1 ratio. Tabulated values 
are only for the Si XII A6. 718 line. 



applying a simple scaling factor to the former such that the 
(local) continua have the same intensity. As the scaling factor 
is wavelength dependent, we divided the two spectra into sev- 
eral bands and applied different factors in each case. The two 
examples shown in Fig. 4 (and all others we have examined) 
illustrate that, except for a luminosity scaling, the two absorp- 
tion line spectra are indistinguishable within the uncertainties 
(^10% for the low-state spectrum and about twice that for the 
high-state spectrum). The somewhat weaker looking lines of 
low ionization species in the high-state spectrum (e.g. around 
7A) are well within the noise. Table 1 lists the EWs measured 
for several key absorption lines (see §3.2.3) in both the high- 
and low-state spectra. Again, the small differences between 
the two spectra are well within the measurement errors. We 
proceed under the assumption that there are no variations in 
the line EWs between the two states. We comment on the 
implications of this in §4. 1 . 

3.2. Spectral modeling 

The second and complementary approach for investigat- 
ing the spectral changes is by detailed modeling of the ab- 
sorbing and emitting gas. The underlying idea is that the 
gas is photoionized by the central X-ray source, and that the 
observed spectra represent its physical state during the two 
states. The principles and the ingredients of such modeling 
were outlined in Netzer (1996) and previous applications to 
the case of NGC 3783 were discussed by Kaspi et al. (2001). 
A very recent analysis is provided by K03 who describe a 
detailed model composed of two absorbers that are different 
from those suggested by Kaspi et al. (2001). Below we com- 
ment on the similarity and differences between our work and 
the new K03 paper Here we summarize our basic method and 
explain its application to the high and the low-state spectra of 
NGC 3783. 

3.2.1. General method and model ingredients 

The X-ray gas is assumed to be photoionized by a central 
source, and in photoionization and thermal equilibrium (see 
§4 for discussion of the last point). Modeling was performed 
using ION2003, the 2003 version of the code ION (Netzer 
1996). The code includes all relevant atomic processes, and 
computes the ionization and thermal structure of the gas along 
with the intensities and EWs for more than 2000 X-ray lines. 
The code is able to consider various geometries, from a sin- 
gle cloud to a multi-component, expanding atmosphere. The 
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Fig. 4. — Left: High-state (red) vs. low-state (blue) 6-7 A spectra of NGC 3783 showing the great similarity in absorption line EWs (cf. table 1). The low-state 
spectrum was multiplied by 1.5 to match the flux of the high-state continuum. For clarity, error bars are only plotted for the low state spectrum. Right: Same for 
the 15.6-18 A range except for the different scaling of the low-state continuum. 



Table 2. Assumed composition 



Element 


Relative abundance 


H 


1.0 


He 


0.1 


C 


3.7 X 10-"* 


N 


1.1 X 10-"* 





5.0 X K)-'* 


Ne 


1.0 X ur"* 


Mg 


3.7 X 10"' 


Al 


3.0 X 10-*^ 


Si 


3.5 X 10-5 


S 


1.6 X 10"' 


Ar 


3.3 X 10"*' 


Ca 


2.3 X lO-*" 


Fe 


4.0 X 10-5 



basic parameters of the model are the gas density (assumed 
to be in the range 10^-^ cm^^ - see justification in §4), the 
hydrogen column density (A^, in units of cm^^), the covering 
factor, the gas composition, turbulent motion (§3.2.2) and the 
oxygen ionization parameter (Uox) defined over the energy 
range 0.54-10 keV. As explained by Netzer (1996), and dis- 
cussed further in George et al. (1998), this choice of ioniza- 
tion parameter (compared with, for example, Ux defined over 
the 0. 1-10 keV range) gives the most meaningful definition of 
the ionization field of X-ray photons for AGN. Regarding the 
covering factor, we distinguish between the emission covering 
factor (Q./4n) applicable to the emission line gas, and the ab- 
sorption (line of sight) covering factor which can be different. 
The "solar composition" used throughout this work is given in 
Table 2 (note the reduced oxygen abundance compared with 
older estimates). 

The incident continuum is taken to be the broken power- 
law defined in Kaspi et al. (2001, Table 4). The only changes 
we have experimented with are related to the slope of the 0.1- 



50 keV continuum. As discussed below, the UV part of the 
spectrum can be different from the one assumed here with 
important implications for the UV absorption lines. This will 
be investigated in a forthcoming paper The effects of the UV 
continuum on the strongest features seen in the X-ray band 
are of far less importance. Thus we consider this SED to be 
adequate for the present analysis. 

The models calculated here are entirely self-consistent and 
are not simple attempts to fit the observed spectra by measur- 
ing line EWs and deducing column densities for the differ- 
ent ions. We have searched for the combination of physical 
components that can be produced in nature in an environment 
where low density gas is exposed to a typical AGN contin- 
uum. These components were then combined in a realistic 
manner, taking into account screening, attenuation of the ra- 
diation field, etc. 

3.2.2. Multi-component models 

We have examined the hypothesis that the spectral changes 
observed arise purely as a result of variations in the opacity 
of the absorbing gas which are caused solely by changes in 
the intensity of the ionizing continuum. Thus, for this experi- 
ment, we assume the shape of the SEDs is the same for both 
the low- and high-states and only differ in total luminosity. 
We start by calculating a variety of models in order to mimic 
the low-state spectrum. Each model is made of several emis- 
sion and absorption components. The absorbers are specified 
by Uox, N and the absorption (i.e. line of sight) covering fac- 
tor The latter is assumed to be the same for all absorbers (but 
see also §4). The absorbers are assumed to be aligned such 
that the observed spectrum is the result of the intrinsic contin- 
uum passing through them all. The emission components are 
specified by their Uox and (that are not necessarily the same 
as those of the absorbers), their covering factor and whether 
or not they are occulted by the absorbing gas. 

The dynamics and kinematics of the absorbing gas are im- 
portant factors in comparing the data with the model. Fol- 
lowing numerous UV observations, and our analysis in Pa- 
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per I, we assumed internal motion in the gas which is much 
larger than the thermal motion. This "turbulent velocity" is 
assumed to have a Doppler profile and the velocity quoted 
is the Doppler b parameter. Paper I showed that all absorp- 
tion lines with good S/N and sufficient resolution (the lines 
of O VII, Ne X, Mg XII and Si XIV) can be characterized by 
two kinematic components. The central velocities of these are 
between -400 and -600 km s"^ and between -1000 and -1300 
km s^^ relative to the systemic velocity. A good represen- 
tation of the observations can be obtained with two Doppler 
profiles each with b 250 km s~^. 

The relative EW and covering factor of the two velocity 
components are critical to our modeling of the source. The 
observations show that the EW ratio of the velocity compo- 
nents is about 1:0.7 (lower:higher) in all lines with sufficient 
S/N (see below). Since this ratio is seen in several saturated 
lines, it can be interpreted as the result of different absorp- 
tion covering factors for each component. In this case, the 
lower- velocity component has a covering factor of 0.8-1.0 
and the larger-velocity component a covering factor of 0.6- 
0.8. However, the composite profiles also include unsaturated 
lines that seem to have similar shapes. In this case, the dif- 
ferent EWs are due to different column densities. We cannot 
distinguish between the two possibihties since we do not have 
high-quality, high-resolution profiles for many weak lines. 
Given these uncertainties, we investigated two cases where 
the relative column densities for the two velocity components 
are in the ratio of 1:0.7. In one case the covering factor of all 
absorbers is unity and in the other case it is of order 0.8. 

In summary, each ionization component of the models pre- 
sented in this paper is made of two kinematic components. 
For the outflow velocities we chose -500 and -1000 km s^\ 
and assumed a turbulent velocity of 250 km s"' for all com- 
ponents. This implies that the total derived column density 
of a certain ionization component is the sum of the colunon 
densities in the two kinematic components, even for saturated 
lines (since there is very little velocity overlap between the 
two). For brevity, in the rest of this paper we just quote the 
total column densities. Given these assumptions, we used the 
EWs measured in Paper I, combined with a few new mea- 
surements, to obtain the various column densities and optical 
depths. The more important lines that were used to constrain 
our models, and their adopted column densities, are listed in 
Table 3. 

3.2.3. The silicon and sulphur line method 

A major clue for the conditions in the absorbing gas is ob- 
tained from EW measurements of various lines in the 5-7.1 A 
band. This wavelength range contains lines from Si VII to 
Si XIV as well as the strongest Unes of S XV and S XVI. Most 
of these lines are not blended. The range of ionization and ex- 
citation is very large and represents the ionization of almost 
all line-producing ions in the Chandra spectrum. The column 
densities deduced from these EWs, assuming the two compo- 
nent profiles, are given in Tables 1 & 3. The atomic data for 
the lines are known either from standard calculations (f- values 
for the H-like and He-like transitions) or from the recent work 
of Behar and Netzer (2002) discussing the inner-shell lines of 
silicon, sulphur and other elements. 

We have developed a simple algorithm to compare the mea- 
sured optical depths of all the observed silicon and sulphur 
lines with the results of the photoionization calculations. We 
first compute a large grid of models with calculated optical 
depths for a large range of ionization parameters and colunon 
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"Log of column density in units of cm ^. 

''Column densities and uncertainties derived from the EWs in Paper 1. The 
EW was divided in the ratio 1 :0.7, then the column density corresponding to 
the 1/1.7 part was computed assuming a Doppler width of 250 km s^'. The 
derived column density was multiplied by 1.7 to result with the total column 
density listed here. Uncertainties were calculated from the EWs uncertanties. 
See text for details. 

"^Saturated line, lower limit only. 



densities. We then pick up to four models at a time and com- 
pare the combined optical depth for the chosen set of lines 
with the values deduced from the observations for the larger 
column density (lower outflow velocity) kinematic compo- 
nent. The best combination of models is obtained by mini- 
mizing the differences between the observed and the calcu- 
lated optical depths in all 10 lines (the 9 silicon and sulphur 
lines listed in Table 1 plus a line of Si VII). The result is a list 
of up to four models, with various column density and Uox, 
whose combination gives the best match to the observed col- 
umn densities. Using this procedure we find that at least three 
ionization components (i.e. six kinematic components) are re- 
quired to fit the data to within the observational accuracy. The 
reason is the very large difference in ionization between the 
lowest (e.g. Si VII and Si VIII) and the highest (e.g. S XVI) 
ionization lines. This requires at least one component to fit the 
EWs of the lowest ionization fines, one to fit the intermediate 
ionization lines and one for the highest ionization lines. 

Experimenting with various modifications of the method 
suggests that there are several combinations of three or four 
ionization components that give similar quality fits to the ob- 
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Fig. 5. — The three generic components that were used to model the low- 
state spectrum of NGC 3783. All components are shown on the same nor- 
malized scale where 1 is the incident continuum level. 



Fig. 6. — A blown up version of Fig. 5 showing the 4.5-7 A region used 
to constrain the model parameters (see text). Note that all three models are 
required to explain the large range of ionization (from Si VII to S XVI.) 



served optical depths of the silicon and sulphur lines. All 
combinations share the same general properties — they all 
require at least three very different sets of physical condi- 
tions (i.e. ionization parameters) and three similar column 
densities. A generic model with the required properties is 
the following three-component model (all column densities 
refer to the total column of the two kinematic components): 
a low ionization component with log{Uox) = —2.4 ±0.1 and 
log(A^) = 2L9±0.1, a medium ionization component with 
\ogiUox) = -1.2 ±0.2 and log(A^) = 22.0 ±0.15, and a 
high ionization component with log{Uox) = —0.6 ±0.2 and 
log(A^) = 22.3 ± 0.2. The uncertainties on Uqx and in the 
case of the low ionization component arise primarily from the 
uncertainty on the slope of the ionizing continuum. Relax- 
ing this constraint allows a somewhat smaller column density 
(see §3.2.4). The larger uncertainties on the parameters of the 
higher ionization components are due to the fact that gas with 
such properties produce similar EWs for medium (e.g. Si XI) 
and high (e.g. Si XIV) ionization lines over a larger range of 
ionization parameters and column densities. The three the- 
oretical ionization components are plotted over a large wave- 
length range in Fig. 5, and a more detailed view of the 4.5-7 A 
range is shown in Fig. 6. 

3.2.4. Intrinsic continuum and global model properties 

The next step is to test the combination of components that 
fit the optical depths of the 5-7.1 A silicon and sulphur ions 
best over other wavelength bands. This requires the shape of 
the underlying ionizing continuum to be defined more pre- 
cisely. The slope of the short wavelength continuum can 
be measured directly if bound-free absorption is negligible. 
However, Fig. 5 clearly shows that two of the generic ab- 
sorbers are characterized by large opacity at wavelengths as 
short as 4-5 A. Thus a unique determination of the 2-5 A con- 
tinuum slope (the part available for the grating observations) 
is not trivial. More specifically, the exact O VII column den- 
sity in the lowest-ionization component strongly influences 
the short wavelength continuum shape, and the uncertainty on 



this column gives rises to an uncertainty on the slope. 

Given the above constraints on the model properties, and 
the observed 2-5 A continuum, we have experimented with 
various power-law continua with the requirement that both 
the short wavelength continuum and the measured EWs of 
the silicon and sulphur lines, are within the observational un- 
certainties. All fits were performed on the low-state spec- 
trum. We found that a single power-law (where the number 
of photons per unit energy, E, is proportional to E^^) with 
F = 1.65 ±0.15 is a good overall representation of the low- 
state spectrum. We note that the limited data quality cannot 
exclude the possibility of some steepening of the intrinsic con- 
tinuum at long {X > 20A) wavelengths. 

We note in passing the different continuum slope adopted 
by K03 in their recent modeling of the source. These authors 
did not consider the two spectral states found here and fitted 
the full Chandra data set. As a result, the SED they adopted 
is based on the combined (high- and low- state) spectrum, and 
hence is different from ours. This is also the reason why their 
powerlaw slope is flatter than the one adopted by Kaspi et al. 
(2001), who fitted the high state only. K03 also include in 
their SED a low-energy blackbody component. 

Another uncertainty on the intrinsic SED arises from the 
fact that some of the high energy photons are likely to be 
produced far away from the source, due to Compton scatter- 
ing and reflection of the central continuum radiation. This 
can be very noticeable in faint sources, like NGC 3516 dur- 
ing 2000 (see Netzer et al. 2002). Scattering by gas in a 
low state of ionization is strongly wavelength dependent, and 
for NGC 3783 can amount to a change of ~ 0.1 in the de- 
rived value of r. As a result, we might have overestimated 
the column density of the lowest ionization component (the 
only component with significant influence on the slope of the 
short wavelength continuum) which would results in overesti- 
mating the column densities of low ionization species such as 
O V and O VI. §4. 1 contains some discussion about the low- 
ering of the column density in our low Uox component and a 
more detailed examination of this effect is deferred to a sepa- 
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rate paper (George et al. 2003, in preparation). 

Finally, a broad, relativistic iron Ka line can also be impor- 
tant in affecting the short wavelength continuum slope. The 
presence or absence of such a feature is an open question. 
Kaspi et al. (2001) see no evidence for this feature in the 
2000 Chandra observation. Given the uncertainties discussed 
above, and the large amount of absorption at almost all wave- 
lengths, we cannot resolve this issue even with our much im- 
proved data. 

Given the chosen continuum and the three ionization com- 
ponents, we have calculated a combined multi-component 
model that covers the entire wavelength range. As stated 
above, the model includes two kinematic components for ev- 
ery ionization component. The theoretical spectrum is com- 
pared with the observed low-state spectrum in Fig. 7. The re- 
sults are very good. In fact, we could not find another combi- 
nation of models (three or more absorbers) that is significantly 
superior. Not only does the model provide a good overall fit 
to many lines of all elements, but it also predicts the strengths 
of the deep bound-free edges correctly. Thus, the constraints 
on the 5-7.1 A silicon and sulphur lines are enough to com- 
pletely specify a model for the entire spectrum. 

3.2.5. Deficiencies in the Modeling 

As demonstrated above, the low-state absorption spectrum 
of NGC 3783 can be satisfactorily explained by our combina- 
tion of line-of-sight absorbers. The agreement with the ob- 
servations is very good for all lines with A < 15 A. However, 
there is a noticeable discrepancy at around 16A. At this wave- 
length range there are two main sources of opacity - the O vn 
bound-free continuum (edge at 16.7 A), and the iron M-shell 
UTA feature (Behar, Sako & Kahn 2001). Our photoioniza- 
tion calculations suggest a disagreement between the two in 
the sense that the strongest iron lines observed are due to ions 
indicative of a lower level of ionization than that indicated by 
the strongest oxygen lines observed. Specifically, the peak 
UTA absorption corresponds to Fe VIII-Fe X. For such a level 
of ionization, our models predict most of the oxygen being 
in the form of O III-O VI. However, the observations indicate 
that most of the oxygen bound-free absorption is due to O VII. 
The possible origin of this discrepancy is discussed in §4. 

Regarding the spectrum beyond 19 A, the comparison be- 
tween the model and the observations is limited by the poorer 
S/N. While the model shown in Fig. 7 is consistent with the 
data, the predicted depths of several absorption Unes in this 
range (in particular those due to O V and O VI) seem to be 
stronger than actually observed. This may be a real short- 
coming of the model (see comments in §3.2.4 and in §4). 
It may also be the result of the assumed absorption cover- 
ing factor. As already mentioned, line-profile analysis sug- 
gest a mean (over the line profile) covering factor of 0.7-1.0, 
and all models considered so far assumed the extreme case 
of an absorption covering factor of unity. The "leakage" of 
the incident continuum makes httle difference at short wave- 
lengths except for a need to somewhat decrease the assumed 
F. However, leakage can influence the comparison at long 
wavelengths much more. This is illustrated in Fig. 8 which 
shows a model with F = 1 .5 and a line-of-sight covering factor 
of 0.85. Indeed, the agreement at long wavelengths is much 
better. The differences between the two assumed slopes and 
covering factors are well within the model uncertainties. 

We have attempted to model the high-state spectrum by as- 
suming the same ionization components and by changing the 



luminosity and the ionization parameter by the factor inferred 
from the observed increased luminosity in the short wave- 
length continuum (~ 1.5). The result is a very poor fit to 
the long wavelength continuum luminosity and to the inten- 
sities of many fines. The reason is that such a small lumi- 
nosity change results in an opacity change which is too small 
to account for the very large difference (factor ^ 4) between 
the low and the high state spectra at long wavelengths. This 
means that the high and the low state spectra are inconsistent 
with the assumption of a simple response of the gas to contin- 
uum luminosity variations. 

There is a simple and satisfactory solution for the high-state 
spectrum that involves the appearance of an additional con- 
tinuum component. This second, long wavelength component 
(sometimes referred to as a "soft excess") appears only during 
high state. We have therefore modeled the high-state contin- 
uum by two powerlaws: the low-state power law (F = 1 .6) and 
a steeper component with F = 3.2 which dominates at long 
wavelengths. We find a satisfactory solution for the high-state 
spectrum when the F = 3.2 component emits 40% of the flux 
of file F = 1 .6 component at 12.398 A (1 keV). (We note that 
this is not a unique combination, and fits of similar quality 
can be obtained for other combinations of slope and normal- 
ization.) The two continua are shown in Fig. 3, along with the 
high and the low state spectra. Given the two component con- 
tinuum, we find a very good agreement for the high-state line 
and continuum spectrum by assuming no opacity variations. 
The quality of the fit is similar to the low-state fits shown in 
Figs. 7 & 8. 

To further iUustrate this point in Fig. 9 we show the high- 
state spectrum divided by the low- state spectrum (binned to 
reduce the noise) along with the ratio of the two assumed con- 
tinua. The diagram shows that any remaining spectral features 
are entirely consistent with the noise. The S/N in this diagram 
is not high enough to completely rule out some opacity-Uke 
variations at long wavelengths. In particular, there is some ex- 
cess emission near the O VII and O VIII recombination edges 
that may hint to extra emission in the high-state spectrum 
(note that we do not expect the emission features to disap- 
pear by this division). However, such opacity variations must 
contribute very Uttle to the spectral variations observed at long 
wavelengths. 

Finally, the emission-line spectrum of NGC 3783 can be ex- 
plained by assuming X-ray emitters with the same properties 
found for the ionized absorbers. The emitted line photons are 
probably observed through the absorber, as indicated by the 
combined emission-absorption profiles of several resonance 
lines. The model is problematic at around 2 1.8 A, where it 
fails to reproduce the O VII intercombination emission line. 
The emission covering factor required is about 0.1 (see §4.1) 
consistent with the value of approximately 0.05-0. 15 obtained 
by Behar et al. (2003). 

4. DISCUSSION 

This section discusses the main results of the new analysis 
of all the Chandra observations of NGC 3783 performed to 
date. The reader is referred also to an appendix where we 
give a detailed comparison with the recent work of Krongold 
et al. (2003). 

The main new findings of our work are the large increase 
in flux between the low- and high-states, which is not accom- 
panied by any significant changes in the properties of the ab- 
sorption. The central source in NGC 3783 fluctuates between 
two such flux states on a 20-120 days time scale, and the best 
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Fig. 7. — A compai'ison of the low state spectium with the six component model (three ionization components each split into two kinematic components). Note 
that most of the apparent discrepancy at around 5-6 A is due to a known calibration problem (H. Marshall, private communication). 



description of the observed variations is the appearance and 
disappearance of a "soft excess", low energy continuum com- 
ponent. These lead to important conclusions regarding the 
nature of the continuum source. Other important conclusions 
are the multi-component nature of the absorber and the large 



range of ionization. Below we discuss some interpretations. 

4.1. Physical properties of the absorbing gas 

We found that three ionization components with different 
properties and a large range of ionization, are required to fit 
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Fig. 8. — Model vs. observation for the long wavelength spectrum assuming F = 1.5 and a line-of-sight covering factor of 0.85 for all absorbers. 
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Fig. 9. — A ratio spectrum of NGC3783: Points with error bars are the 
result of dividing the high-state spectrum by the low-state spectnim. The 
binning is 0.04A for A < UA and 0.1 A at longer wavelengths. The smooth 
line is the same division applied to the two continua shown in Fig. 3. 



the spectrum of NGC3783. Each of those components is 
split into two kinematic components, as implied by the ob- 
served line profiles reported in Paper 1 . The three "generic ab- 
sorbers" of our model are listed in Table 4 and shown in Figs. 
5 & 6. As explained above, the number of real components 
and the exact values of Uqx and are not very meaningful 
since several somewhat different combinations of ionization 
components give equally good fits. Rather, the values chosen 
here are representative, "average" ionization parameters, and 
the total column densities required. Our modeling of the low- 



Table 4. Parameters of the Three Generic Models 



Model 


\og{Uox) 


log(/V) 


Low ionization 


-2.4 


21.9 


Intermidiate ionization 


-1.2 


22.0 


High ionization 


-0.6 


22.3 



State spectrum assumes the underlying X-ray continuum is a 
single power-law. However, the S/N of the long wavelength 
spectrum is not high enough to completely rule out the pres- 
ence of an additional softer continuum component, similar to 
the one seen during the high states but of much lower lumi- 
nosity. We also note that there is little (if any) response of 
the absorbing gas to the short (~ 2 days) time scale contin- 
uum variations. Thus our measurements are models represent 
some time-averaged spectral properties. 

The main difference between the present model and the 
one given by Kaspi et al. (2001) is the presence of the low- 
est ionization component. Only the first Chandra observa- 
tion (obtained in 2000 January) was available to Kaspi et al., 
and hence was insufficient to reveal the two spectral states. 
Furthermore, this relatively short observation resulted in the 
S/N at A > 18A which was insufficient to reveal the true 
shape of the continuum. Hence Kaspi et al. naturally (but 
incorrectly) assumed the underlying continuum was a single 
power-law for this high-state spectrum. As a result, their 
preferred model contained a component with a strong O VII 
bound-free edge, but very little absorption at longer wave- 
lengths (i.e. at A > ISA the underlying continuum was "re- 
covered"). Our new observations do not show such a recovery. 
In fact, the absorption in the long wavelength region is dom- 
inated by features of K-shell carbon and nitrogen and L-shell 
nitrogen and oxygen. 

The presence of a lower ionization component is also a key 
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point in the K03 work. In their study they assume two ab- 
sorbers with log{Uox) = —2.77 (i.e. a factor of about two 
lower than our lowest ionization component) and log{Uox) = 
— 1.23 (very similar to our imtermediate ionization compo- 
nent). Their fit does not require a third, very ionized com- 
ponent. Our reason for suggesting this higest-ionization 
(log{Uox) = —0.6) component is based on the fit to the lines 
of Si XIV, S XV and S XVI. When this component is not in- 
cluded, the discrepancy in column densities of these ions is 
a factor 2-3. Several of the lines of these ions are in fact 
missing from the K03 analysis. Moreover, some of the high- 
est -ionization lines are badly saturated, (e.g. Si XIVA6.18A). 
Thus a comparison of the observed and the calculated EWs of 
those lines, as performed by K03, is not sufficient for evaluat- 
ing the ionic column densities. 

Perhaps the largest observed discrepancy between model 
and observation is due to absorption by the iron M-shell UTA. 
This feature has been observed in several other AGN (Steen- 
brugge et al. 2003 and references therein) with similar shape 
and central wavelength. Its smooth shape is due to the large 
number of absorption lines of many iron ions and the wave- 
length of largest absorption is ionization dependent. Our mod- 
eling of this spectral region includes all the lines in the Behar 
et al. (2001) calculations (i.e. not only the shorter published 
list of lines). The central observed wavelength of the fea- 
ture in NGC 3783 is at around 16. 4A suggesting that most 
of the contribution is due to Fe VIII-Fe X. For the assumed 
SED, this corresponds to gas whose dominant oxygen ions 
are O lU-O VI. In our modeling, most long wavelength ab- 
sorption is due to the component with \og{Uox) = —2.4 (Fig. 
5). The dominant ions in this components are O VI-O VII and 
Fe X-Fe XII. This corresponds to peak UTA opacity at around 
16.1 A. As shown in Fig. 7, the model fits the observed oxy- 
gen lines and continua reasonably well but badly misses the 
position of the UTA absorption. 

We have experimented with various other models to try to 
eliminate this discrepancy. For example, an overabundance 
of iron relative to oxygen may decrease the level of ioniza- 
tion of iron due to the increased iron opacity. Experiment- 
ing with iron abundance which is three times larger resulted 
with changes that are not large enough to explain the observed 
discrepancy. We have also experimented with models that in- 
clude several low ionization components, instead of the single 
low Uox component shown here. This, again, gave very lit- 
tle improvement. It seems that the apparent conflict between 
the oxygen and iron ionization cannot be resolved by these 
models and we suggest two alternative explanations. 

1. The absorbing gas in NGC 3783 may not be in ioniza- 
tion equilibrium due to its low density and the rapid 
flux variations on time scales that are shorter than the 
recombination and ionization time scales. This is prob- 
ably not very important for the short wavelength flux 
since its variability time scale is short and the ampli- 
tude not very large. The gas responding to this con- 
tinuum is probably at some mean level of ionization. 
This however is not the case for the soft excess con- 
tinuum which varies on a much longer time scale with 
a much larger ampUtude. This can contribute, signifi- 
cantly, to the ionization of the lower ionization species. 
Different ions react on different time scales and the gas 
may never reach an equilibrium. Since all our mod- 
els assume steady state gas, they may not be adequate 
to describe the absorbing gas properties. In particular. 



the iron and oxygen recombination times may be dif- 
ferent enough to result in gas where iron is less ionized 
or oxygen is more ionized compared to the equilibrium 
situation. The complicated issue of time dependent ion- 
ization is beyond the scope of this paper. 

2. The ionization balance of iron depends critically on low 
temperature dielectronic recombination (DR) rates that 
are not well known for the iron M-shell ions. Experi- 
ence with L-shell iron ions, whose DR rates have been 
measured and computed recently (Savin et al. 2002 and 
references therein), shows that in almost all cases the 
previous low temperature DR rates were consistently 
smaller than the newly calculated values. Assuming a 
similar effect in the M-shell ions, we can envisage a sit- 
uation where, in the absence of realistic DR rate, the en- 
tire ionization balance of iron is shifted toward higher 
ionization. Thus, more realistic low temperature DR 
rates may bring models to better agreement with the ob- 
servations. We note that low temperature DR rates are 
also not available for the magnesium and siUcon ions 
that we have modeled. For some of these ions (Mg VIII 
and Mg IX) these are probably not very important and 
for others (the relevant silicon ions) the situation is less 
clear (D. Savin, private communication but see also Gu 
2003). Nevertheless, since we do obtain good agree- 
ment between the calculated EWs of these lines and the 
EWs of other lines observed in the spectrum, we did not 
consider changing those rates in the calculations. 

It is interesting to note that the K03 model assumes a lower 
level of ionization for all elements, including iron, and hence 
provides a better fit to the UTA. However, their lower ion- 
ization component fails to fit the Si X and Si XI lines around 
6.8A. Our model was designed to give a very good fit over the 
5-7A wavelength range and hence produces a lower quality 
fit to the UTA feature. Thus, there seem to be no satisfac- 
tory solution based on a single low-ionization component that 
explains all these features. 

An interesting and somewhat problematic issue is the abun- 
dance of O VI, which can also be clearly observed in the 
far UV (Gabel et al. 2003a; Gabel et al., 2003b). Paper I 
found even the Ka line of O VI undetectable with the cur- 
rent S/N. However, revisiting this measurement more care- 
fully, and with more accurate atomic data, we can now detect 
several O VI lines. From two such lines, both of which are 
probably saturated, we have obtained a lower Umit of lO'^^ 
cm~-^ on the column density of O VI. This is much higher than 
the value (lo'^ ''^''-^ cm"^) obtained from the XMM-Newton 
RGS spectrum (Behar et al. 2003). Moreover, the values pre- 
dicted from our model (10^^^^ cm^^) is even larger. Such a 
large value is problematic for several reasons. It may be in 
contradiction with the UV observations of the source (see be- 
low), it is much larger than the values obtained from the RGS 
observations, and it results in a prediction of a strong absorp- 
tion line at around 21.87 A which is not seen in the spectrum 
(Fig. 7, the absorption line underneath the O VII intercombi- 
nation line). 

An explanation of the discrepancy between the values for 
the O VI column density measured by Chandra and XMM- 
Newton is the crowded region of the spectrum where O VI 
Ka resides. In particular, the lower spectral resolution of the 
XMM-Newton RGS makes it extremely difficult to resolve the 
weak O VI absorption line at 22.01 A in the presence of the 
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Fig. 10. — The fractional ionization of oxygen, magnesium and silicon 
ions in a low density gas of our chosen composition exposed to the NGC 3783 
low-state continuum. The curves show the mean ionization of various ions 
averaged over a column density of 10^' cm"^. Note the great similarity in 
the ionization of O 'VI, Mg "VIII and Si "VIII (all plotted with dotted Hues). 



that K03 analyse the combined (high- and low-state) spec- 
trum, which is then interpreted in terms of a different amount 
of bound-free absorption. Another reason is the somewhat 
larger oxygen abundance assumed by K03 (the oxygen col- 
umn densities of the low-ionization components in the two 
papers actually differ only by a factor 1.35). We also note that 
the inner-shell O VI lines that we have used to determine the 
limit on the column density of this ion are missing from the 
K03 model. This makes a noticeable difference in the spectral 
fitting at around 21 A. 

Finally, we have looked for the lowest O VI column density 
allowed by the observations if we remove the constraint im- 
posed by the short wavelength continuum slope (§3.2.4) and 
keep the ionization parameter unchanged. This means find- 
ing the lowest column density for the log([/ox)=-2.4 compo- 
nent which is still consistent with the measured EWs of the 
5-7.1 A silicon and sulphur lines. We found this limit to be 
log(A^) > 21.7 for the low-state spectrum, and log(A^) > 21.6 
for the high-state spectrum. These lower limits on the column 
densities are consistent with the uncertainties given in Table 
1 . Under these assumptions, the column density of O VI in 
the low-state spectrum can be as low as lO'^^^ cm^^. This 
column is large enough to produce saturated O VI lines at 
21.01 and 21.87 A. It also requires a much harder continuum 

(r=i.4). 



bright O VII forbidden emission line (22.10 A). We also wish 
to note that with the lack of laboratory measurements for the 
O VI wavelengths, high-resolution astrophysical spectra such 
as the present one constitute the best determination of these 
wavelengths currently available. The Chandra data, along 
with the assumption that O VI is outflowing at the same ve- 
locity as O VII, places the O VI Ka line at a rest-wavelength 
of 22.01 ± 0.01 A. This value is in agreement with our model 
wavelength of 22.005 A, calculated with the HULLAC atomic 
code, as well as with the value measured in NGC 5548 (Steen- 
brugge et al., 2003), but is somewhat different from the value 
of 22.05 A calculated by Pradhan (2000) using the R-matrix 
method. For more information see Behar and Kahn (2002). 

A O VI column density less than lO'^ cm^^ is also incon- 
sistent with theoretical predictions based on modelling the 
abundances of Mg VIII, Si VIII and O VI. This is illustrated in 
Fig. 10 where we show calculations for the fractional abun- 
dance of oxygen, silicon and magnesium ions as a function of 
Uox assuming an incident continuum similar to the one used 
here. The diagram shows the very similar fractional abun- 
dance of O VI, Mg VIII and Si VIII as a function of ioniza- 
tion parameter. Using the measured EW of the Si VIIlA 6.998 
line, which is an unblended, easy-to-measure feature, we can 
use Fig. 10, and the assumed silicon-to-oxygen abundance 
ratio, to derive an estimate for the column density of O VI. 
This number is very close to the model prediction. The sit- 
uation is more complicated regarding magnesium since sev- 
eral of the relevant lines (around 9.4A) are blended with neon 
lines. However, we could obtain an estimate of the EW of the 
Mg VIII A9. 506, 9.378 lines that are good enough to constrain 
the Mg VIII column density. The O VI column density based 
on this measurement is, again, much larger than lO'^ cm^^ 
and in good agreement with the model prediction. 

A comparison with the K03 work shows a disagreement by 
a factor of two in the total (hydrogen) column density of the 
lowest-ionization absorber. One reason for the difference is 



4.2. Outflow velocity and covering factor 

Some of the parameters our model are based on measured 
profiles of several absorption lines that were shown in Paper I 
to include at least two components with different ouflow ve- 
locities. The measurements in Paper 1 also show that the rela- 
tive EWs of the two components seems to be ionization inde- 
pendent and that most of the lines used for the profile analysis 
were heavily saturated. This suggests that the apparent opti- 
cal depths of those lines are determined by the covering fac- 
tor rather than the line opacity, and hence similar to the case 
seen in many UV absorption systems in AGN (e.g. Barlow, 
Hamann & Sargent 1997;Arav et al. 1999). The covering fac- 
tor of the lower- velocity components of the O VII, Ne X and 
Si XIV lines is in the range 0.8-1 .0. The covering factor of the 
larger-velocity component of the same lines, assuming satu- 
rated profiles, is 0.6-0.8. We also note that the large number 
of iron lines near the UTA center gives an independent esti- 
mate of the covering factor which is at least 0.85. 

Fig. 8 shows that models with a covering factor less than 
unity provide fits to the data of similar quality. Using the sil- 
icon and sulphur lines (as in §3.2.3), we find three ionization 
components are again required, each with parameters similar 
to those listed in Table 4. The only difference is the need 
to assume a somewhat harder ionizing continuum. The fit 
to the long wavelength lines is in fact somewhat better than 
that produced by the full covering models. However given 
the low S/N at those wavelengths, this is not a strong con- 
clusion. Thus, the data-model comparison cannot constrain 
covering factor beyond the actual observations. The most im- 
portant conclusion is that the covering factor is similar in low 
(e.g. O VII and the iron UTA) and in high (e.g. Si XIV, Ne X) 
ionization lines. 

4.3. Density, location and thermal stability of the absorbing 

gas 

Our measurements and analysis are consistent with no vari- 
ations in the absorbers' properties on time scales of 1^ days. 
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There are also indications for no spectral changes over much 
longer periods, perhaps a few months. For instance we find no 
evidence for any significant, narrow (< 1 A) spectral features 
in the high- to low-state spectral ratio shown in Fig. 9. The 
reason we cannot constrain the variability of the absorbers 
on short time-scales are the fast, relatively-small amplitude 
fluctuations in the luminosity of the central source (maximum 
amplitude of about a factor 2). As a result of these variability 
characteristics of the illuminating continuum, the lack of sig- 
nificant spectral variations may be due to the ionization level 
of the gas reflecting the time averaged ionizing luminosity, 
rather than very long ionization or recombination times. This 
conclusion therefore depends primarily on the direct compar- 
ison between the low and the high state spectra, and thus on 
the S/N in the lower-quality spectrum used for the comparison 
(the high state spectrum). 

To further examine this point, we have calculated the theo- 
retical spectrum appropriate for the high state continuum un- 
der the assumption that the absorbing gas responds instantly 
to changes in the ionizing continuum. This would be the case 
if the gas along the line-of-sight gas were of very high den- 
sity. The initial conditions are those assumed for the low-state 
continuum, and the change in the SED as described earlier: a 
two component X-ray continuum, where the long-wavelength 
component is assumed to extrapolate down to 40 eV. The in- 
crease in flux between the low- and high states was assumed 
to be a factor 1.5 at 4 A (see Fig. 3). Fig. 1 1 shows two the- 
oretical spectra: the low-state spectrum (in blue) is the same 
model shown in Fig. 7, the high state spectrum (in red) is the 
result of these new calculations. The two are clearly differ- 
ent. In particular the Si VII, Si VIII and Si IX lines are are 
much weaker for the high-state model. This is due to the 
large increase in flux at longest wavelengths, with the largest 
consequences for the lowest ionization component (Uox has 
increased by a factor close to 2 and Ux by a factor of 6.6). 
The variations predicted by this theoretical exercise are much 
larger than those observed (Table 1 and Fig. 4). This illus- 
trates the gas did not simply react to the differences in the 
illuminating continuum between the low- and high-states. An 
even stronger conclusion is obtained from the fact that at long 
wavelengths the differences observed between the two states 
are much larger than can be explained by simple variations in 
opacity. 

Given the above findings, we can place limits on the elec- 
tron density (n^,) and the distance from the central source (D) 
for each absorption component. We used the temperatures 
derived from the models, average recombination times for 
the dominant oxygen and silicon ions, and assuming there is 
no response to continuum variations on a time scale of 10 
days. We find rie < 5 x lO'* cm^^ and D > 3.2 pc for the 
log{Uox) = —2.4 component, rie < 10^ cm^^ and D > 0.63 
pc for the log{Uox) = — 1.2 component , and rig < 2.5 x 10^ 
cm^^ and D > 0.18 pc for the \og{Uox) — —0.6 compo- 
nent. The corresponding masses are 8.1 x lO^C/, 390C/ and 
64C/ solar masses, respectively, where Cf is the absorption 
(47r) covering factor The mass outflow rate is dominated 
by the lowest ionization component and is approximately 
75C/'£Voutflow/(500 kms^') solar masses per year, where £ 
is the radial filling factor of the flow. The very large num- 
ber suggests a short duration outburst rather than continuous 
ejection. We also note that the covering factor can be small, of 
order 0.1. The above distances are in agreement with the find- 
ings of Behar et al. (2003), who have suggested the outflow is 
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Fig. 1 1. — Low state and high state theoretical spectra assuming the ob- 
served variable NGC 3783 continuum and line of sight gas that responds in- 
stantly to continuum variations. The low state model (in blue) is the one 
described earlier and shown in Fig. 7. The high state spectrum (in red) as- 
sumed the high state SED discussed in the text. Such a continuum results in 
much weaker Si VIII and Si IX lines. 



located at a distance of a few pc and is extended beyond 10 pc 
based on a comparison between the soft X-ray emission and 
absorption characteristics of the XMM-Newton data. 
A significant new aspect of our model is that the product 
X r is similar (within a factor <2) for all three absorp- 
tion components. This raises the interesting possibility that 
all the components are in pressure equilibrium (assuming gas 
pressure dominates the total pressure, i.e. radiation pressure 
and turbulent pressure are not important), and they all occupy 
the same volume of space. To test this idea, we have calcu- 
lated thermal stability curves (log(r) vs. \og{Uox /T) for the 
low-state continuum and our assumed composition (Table 2) 
under various assumptions. We have kept the UV SED un- 
changed and varied r(0.1-50 keV) over a wide range. We 
have also investigated the possibility of a line-of-sight atten- 
uation by the absorbing gas. This can be relevant for the sit- 
uation under study since, even in the simplest geometry, two 
of the absorbers do not have a clear view of the central radi- 
ation source. Screening by the lowest-ionization component 
{\og{Uox) = —2.4) is the most important since this gas mod- 
ifies the transmitted spectrum much more than the other com- 
ponents (see Fig. 5). Three such stability curves are shown in 
Fig. 12, along with the location of the three ionization com- 
ponents. The exact location depends on the details of the as- 
sumed SED, the gas metallicity and opacity. 

The interesting feature of the stability diagram is that all ra- 
diation fields considered here result in extended, almost verti- 
cal parts between T ~ 3 x 10^ K and T ~ 2 x 10^ K. The curve 
for F = 1 .6 has a more extended vertical part since the higher 
mean energy of the radiation field results in a Compton tem- 
perature. Screening by the lowest-ionization component in 
NGC 3783 is indeed important, but does not change the main 
result that, given the uncertainties on the SED and compo- 
sition, all three ionization components lie on the stable part 
of the curve. We also note that, while here we only consid- 
ered gas with a single composition, the details of the stability 
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Fig. 12. — Thermal stability curves for a low density gas exposed to the 
low-state continuum of NGC3783. The two curves marked by F are for a 
bare continuum with a 0.1-50 keV slope as marked. The dotted line is the 
stability curve for a gas exposed to the F = 1.8 continuum seen through the 
log(J/ox) = —2.4 absorber. Thick sections mark the locations of the three 
absorbers considered in this work with the corresponding uncertainties on 
the ionization parameters. All ionization components considered in this work 
are situated on thermally stable parts and they all have roughly the same gas 
pressure. 



curve are abundance dependent. Needless to say, very differ- 
ent curves will be obtained for significantly different SEDs. 

We note that K03 reach similar conclusions regarding the 
pressure equilibrium of their two ionization components. 
However, in our model all components are situated on stable 
branches of the thermal stability curve (Fig. 12) while accord- 
ing to K03, their high-ionization component is situated on the 
unstable branch (their Fig. 16). We suspect that the difference 
is in part due to the different SED and metahcity assumed by 
K03, and in part due to the method they used to calculate the 
curve. 

The existence of an additional stable-region along the curve 
of thermal-equilibrium (intermediate between the cold and the 
hot branches) where "warm" {few x 10^ K) can survive is not 
a new idea. It has been discussed by Marshall et al. (1993) 
in their study of the multi-phase medium in NGC 1068, and 
later by many others. Reynolds & Fabian (1995) discussed 
the stability criteria and noted the narrow range of ionization 
parameter allowed for warm absorbers. Hess et al (1997) dis- 
cussed the various parameters governing the stabiUty curve 
(metallicity and more). Komossa and collaborators (e.g. Ko- 
mossa & Fink, 1997; Komossa & Meerschweichen 2001) cal- 
culated many such curves, for different SEDs and metallici- 
ties, and discussed the location of warm absorbers as well as 
the possible Unk between the NLR and the X-ray gas. Kro- 
lik & Kriss (2001) considered multi -phase, warm-absorber 
winds and noted the various locations of stable, over-heated 
and over-cooled gas. Kinkhabwala et al. (2002) re-visited 
the thermal stabihty issue in NGC 1068 and showed that a 
wide range of Uqx is needed to account for the observed X- 
ray emission lines. Chelouche & Netzer (2003) investigated 
the changes in the stability criteria for clouds with large inter- 
nal line radiation pressure. 



Most of these earlier works focused either on the allowed 
parts of the curve from general stability conditions, or on the 
location of a certain absorber observed in a certain source. 
Our work shows that three very different aborbers in a source 
can all occupy extended, stable parts of the curve where the 
gas pressure is roughly the same. This means there may well 
be more than three absorbers spread over the vertical branch 
of the curves in Fig. 12. The three ionization parameters con- 
sidered here may therefore represent some volume-averaged 
properties of the entire cloud ensemble. The column densi- 
ties are therefore the total column densities of a large num- 
ber of such clouds. Moreover, the most-ionized component 
may provide the confining medium for the two other ioniza- 
tion components. Needless to say, real confined outflowing 
clouds are characterised by a more complicated density and 
pressure structure than with the simplified, constant density 
clouds considered here (e.g. see Chelouche & Netzer 2001). 

The pressure equilibrium between all three absorbers in 
NGC 3783 suggest they occupy the same general location. 
This aUows additional constraints to be placed on their density 
and location. Specifically, we require that the distance of this 
region is at least 3.2 pc (the minimum distance of the lowest 
ionization component). However, this zone cannot be much 
further away because of total size and volume considerations. 
Under the assumed conditions, the highest-ionization param- 
eter component has the lowest density, and hence the largest 
dimension. However, its radial extent cannot exceed about 
500 pc (the approximate radius of the narrow-line region), 
and thus its mean density cannot fall below about 10 cm^^. 
Pressure equilibrium then dictates the density of the lowest- 
ionization component, and a maximum distance of about 25 
pc. This limitation is only imposed on discrete clouds — i.e. 
on the mean properties of the gas. Real absorbers may not be 
"clouds", they can cover a large range of distances, and can 
have properties which depend on location. The detailed in- 
vestigation of such models is beyond the scope of the present 
paper. 

4.4. The emission lines 

There is no indication for flux variation in any of the ob- 
served emission lines. However the measurement uncertainty 
on these line is rather large since the strongest lines are ob- 
served at the long wavelengths where the S/N is poor We 
obtain satisfactory fits for most lines by assuming that the 
emission-line gas has the same column density and ionization 
parameter as the absorbing gas. The covering factors required 
for the emitting gas are 0.2-0.3 for the \Qg{Uox) = —2.4 com- 
ponent, 0.1-0.2 for the log(Uox) — —1-2 component, and 
~ 0.1 for the \og{Uox) = —0.6 component. The large range 
is due to the uncertainties in the continuum placement and in 
estimating the fraction of the O VII and O Vin emission lines 
that are being absorbed by the line-of-sight gas (the model 
shown in Figs. 7 and 8 assumes that all the emitted photons 
are seen through the absorbers). 

A clear shortcoming of the model is the fitting of the in- 
tercombination O vn Une. The intensity predicted for this 
line is too small by a large factor (2-4). A large part of the 
discrepancy is due to the O VI A21.87 A absorption line sit- 
uated close to this emission (Fig. 7). This may be related to 
the O VI problem discussed earher. Alternatively it may be a 
consequence of our assumption that all emission-line photons 
are seen through the absorbers. Indeed, a model that assumes 
no absorption of the emitted photons by the line-of-sight gas 
(not shown here) would give a much better agreement to the 
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O VII absorption complex. This would also required a smaller 
emission covering factor Obviously we cannot determine the 
exact geometry using current observations, and the real case 
may well be intermediate between those two extremes. 

We also note on the difference between the absorption cov- 
ering factor (0.8-1.0) and the emission covering factor (0.1- 
0.3). Given the uncertain geometry and the fact that the first 
is a line-of-sight covering factor and the secong a globel (An) 
property, the two are perhaps consistent with each other. We 
did not investigate aU possibilities and cannot comment in de- 
tail on special cases like biconical flows or a torus geometry 
(e.g. Krolik & Kriss 2001). 

Finally, comparing our method with that of K03 we note 
that those authors do not calculate the emission line flux and 
profile in a self-consistent way. Instaed, they add Gaussian 
shaped emission profiles to their calculated absorption pro- 
files (e.g. their Fig. 6). This can result in significant differ- 
ences regarding the line flux and EW. For example, our default 
approach assumes that all emission lines are seen through the 
absorber. In such a case, the entire "blue" wing of many emis- 
sion lines is absorbed. Clearly this affects the EW of the emis- 
sion lines, but has almost no effect on the EW of the absorp- 
tion features. In contrast, the approach adopted in K03 always 
results in a decrease of the EW of the absorption lines. 

4.5. UV absorbers in NGC 3783 

NGC 3783 contains time-variable absorbers in the UV that 
have been described in various earlier publications (Gabel et 
al. 2003a and references therein). Our analysis clearly indi- 
cates that at least one X-ray component contains gas which is 
of sufficiently low ionization so as to produce strong absorp- 
tion lines of C IV, N V and O VI in the UV. 

The predicted EWs of these UV lines depend on two fac- 
tors that have been given little attention in this work: the ex- 
act shape of the UV continuum, and covering factor appro- 
priate for the UV absorbers. A UV continuum which is softer 
than that assumed here can result in higher levels of ionization 
for carbon, nitrogen and oxygen, with little effect on most X- 
ray lines. This is very important for the C IVA 1549 A and 
N VA1240 A lines, but less so for the UV lines of O VI tiiat 
are directly linked to the O VI column density discussed here. 

As for the covering factors, the size of the UV continuum 
source is likely to be much larger than the size of the central 
X-ray source. Thus the covering factor appropriate in the UV 
can be different to that appropriate in the X-ray. This will re- 
sult in saturated UV absorption lines whose EWs are smaller 
than the ones indicated by the column densities derived here. 
However, our upper limit on the density of the low-ionization 
component, combined with its column density, suggest a very 
large line-of-sight dimension (> lO'^ cm). Assuming lateral 
dimension of the same size or large, we conclude that the 
physical dimension of this component is much larger than the 
expected dimension of any likely UV source (e.g. the surface 
of a thin accretion disk). A combined analysis of the UV and 
X-ray results, taking these such points into consideration, is 
in progress 

UV absorbers are also UV emitters, and the low-ionization 
component can also contribute to the observed UV emission 
lines. Our photoionization calculations show that this con- 
tribution is very small. For example, an emission covering 
factor of 0. 1 , similar to the one deduced from the X-ray emis- 
sion lines, wiU produce O VIA 1035 line with emission EW of 
about lA. 



4.6. The long wavelength component 

Perhaps the most interesting result of this study is the ap- 
pearance and disappearance of the soft excess component. 
This broad-band continuum source was seen in two of the 
six observations. On average, these are the two observations 
with the highest luminosity. An increase in softness ratio with 
increasing source luminosity is well known and well docu- 
mented in a number of AGN (e.g. Magdziarz et al., 1998; 
Chiang et al., 2000; Markowitz & Edelson 2001). However 
we are not aware of a softness ratio increase that is not corre- 
lated with a short wavelength flux increase. NGC 3783 seems 
to be the first AGN to show this phenomenon, but we suspect 
that careful spectroscopic monitoring will reveal the same or 
a similar behavior in other sources. This has important conse- 
quences to the continuum production mechanism, as well as 
for the modeling of the warm X-ray gas around the center. 

We do not know the origin of the soft continuum source. 
In particular, we do not have information on its flux at wave- 
lengths A > 30A (the source was not in a high-state in any 
of the published XMM-Newton observations). It may be re- 
lated to a broad-band phenomenon (e.g. due a flaring accre- 
tion disk), or a component covering a narrower band (e.g. a 
single temperature black body). Simple global energy consid- 
eration show that this phenomenon is unrelated to the appear- 
ance of broad emission features like those claimed to be seen 
in at least two sources and interpreted as due to relativistic 
disk lines (Branduardi et al 2001; Mason et al. 2002; but see 
also Lee et al. 2002). 

Finally, we must comment on the possibility of a more com- 
plex behavior. Our analysis is based on the fitting of a single 
powerlaw continuum at low-state and the addition of a soft 
X-ray component during the high-state. However, an equiv- 
alent analysis could be carried out in the reverse order - i.e. 
starting from a pure powerlaw for the high-state and subtract- 
ing a continuum component to explain the low-state spectrum. 
However, as argued earlier, the observed spectral changes at 
long wavelengths are too large to be explained by pure opacity 
variations and it is not at aU clear what other mechanism could 
explain an AGN continuum with flux deficit at long wave- 
lengths. 

5. CONCLUSIONS 

Our detailed measurements and analysis of the 900 ks data 
set of NGC 3783 lead to the following results: 

1. The source fluctuates in luminosity, by a factor ~ 1.5, 
during individual 170 ks observations. The fluctuations 
are not associated with significant spectral variations. 

2. On time scales of 20-120 days, the source exhibits two 
very different spectral shapes denoted here as the high- 
and low-states. The two are associated with different 
softness ratios that seem unrelated to the total X-ray lu- 
minosity. The observed changes in the underlying con- 
tinuum can be described as due to the appearance (in 
the high-state) and disappearance (in the low-state) of 
a soft excess component. The origin of this continuum 
component is not clear. To the best of our knowledge, 
NGC 3783 is the first AGN to show such a behavior. 

3. The appearance of the soft continuum component can 
explain all spectral variations observed within the mea- 
surement uncertainties. There is no need to invoke 
opacity changes between the high- and low-states. This 
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conclusion depends mostly on the S/N in the high state 
spectram. 

4. A combination of three ionization components, each 
split into two kinematic components, provides a good 
description of the intensity of almost all the absorption 
lines and bound-free edges observed. The components 
span a large range of ionization, and have a total col- 
umn of about 4x10^^ cm^'^. The only real discrep- 
ancy between the observed and theoretical spectra are 
for the iron M-shell UTA feature at 16-16.5A. This is 
most likely due to inadequate dielectronic recombina- 
tion rates currently available to use in our calculations. 
The largest other uncertainty is in the column density 
of O VI. 



5. The three generic absorbers discussed in this work have 
very similar values of tig x T . We speculate that the 
absorbers may be in pressure equilibrium with each 
other, occupying the same volume in the nucleus. This 
is the first confirmation of the location of several X- 
ray absorbers on the vertical part of the log(r) vs. 
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6. We have obtained thee lower limits on the gas distance 
from the center, corresponding to our three generic ab- 
sorbers. The limits are 3.2 pc, 0.6 pc and about 0.2 
for the low ionization, intermediate ionization and high 
ionization absorbers, respectively. The pressure equi- 
librium assumption implies distances in the range 3- 
25 pc. 
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